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1—5 Abstract 

Models of tumor growth, now commonly used, present several levels of complexity, 
both in terms of the biomedical ingredients and the mathematical description. The 
simplest ones contain competition for space using purely fluid mechanical concepts. 
^ Another possible ingredient is the supply of nutrients through vasculature. The models 

can describe the tissue either at the level of cell densities, or at the scale of the solid 
tumor, in this latter case by means of a free boundary problem. 

Our first goal here is to formulate a free boundary model of Hele-Shaw type, a vari- 
ant including growth terms, starting from the description at the cell level and passing 
to a certain limit. A detailed mathematical analysis of this purely mechanical model 
is performed. Indeed, we are able to prove strong convergence in passing to the limit, 
with various uniform gradient estimates; we also prove uniqueness for the asymptotic 
Hele-Shaw type problem. The main tools are nonlinear regularizing effects for certain 
porous medium type equations, regularization techniques a la Steklov, and a Hilbert 
duality method for uniqueness. At variance with the classical Hele-Shaw problem, here 
q the geometric motion governed by the pressure is not sufficient to completely describe 

the dynamics. A complete description requires the equation on the cell number density. 

Using this theory as a basis, we go on to consider the more complex model including 
nutrients. We obtain the equation for the limit of the coupled system; the method relies 
^ on some BV bounds and space/time a priori estimates. Here, new technical difficulties 

appear, and they reduce the generality of the results in terms of the initial data. Finally, 
we prove uniqueness for the system, a main mathematical difficulty. 
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1 Motivation and tumor growth models 



In the understanding of cancer development, mathematical modeling and numerical simula- 
tions have nowadays complemented experimental and clinical observations. The field is now 
mature; books and surveys are available, as for example [31 IH 24, 36J. A first class of models, 
initiated in the 70's by Greenspan [31], has considered that cancerous cells multiplication is 
limited by nutrients (glucosis, oxygen) brought by blood vessels. Models of this class rely 
on two kinds of descriptions; either they describe the dynamics of cell population density 
[12] or they consider the 'geometric' motion of the tumor through a free boundary problem; 
see [T5| [lUl |2"o] and the references therein. This stage lasts until the tumor reaches the 
size of ~ 1mm; then lack of nutrients leads to cell necrosis which triggers neovasculatures 
development [15] that supply the tumor with enough nourishment. This has motivated a 
new generation of models where growth is limited by the competition for space [TT] , turning 
the modeling effort towards mechanical concepts, considering tissues as multiphasic fluids 
(the phases could be intersticial water, healthy and tumor cells, extra-cellular matrix . . . ) 
[T3] dH [37] . This point of view is now sustained by experimental evidence [3B] • The term 
'homeostatic pressure', coined recently, denotes the lower pressure that prevents cell multi- 
plication by contact inhibition. 

The aim of this paper is to explain how asymptotic analysis can link the two main ap- 
proaches, cell density models and free boundary models, in the context of fluid mechanics. 
We depart from the simplest cell population density model, proposed in [13], in which the 
cell population density g(x, t) evolves under pressure forces and cell multiplication according 
to the equation 



where p is the pressure field. Pressure-limited growth is described by the term $(p), which 
typically satisfies 



for some pu > (the homeostatic pressure). The pressure is assumed to be a given increasing 
function of the density. A representative example is 



with parameter m > 1. In the free boundary problem to be discussed later, the value q c 
represents the maximum packing density of cells, as discussed in [4"T] . 

Nutrients consumed by the tumor cells and brought by the capillary blood network are a 
usual additional ingredient to the modeling. The situation is described in that case by the 
system 



(1.1) 



d t Q - div(gVp) = q $(p), 



(1.2) 



$'(p) < and $(pm) = 



(1.3) 




(1.4) 



d t g - diy (qVp) = Q$(p,c), 
< d t c — Ac = — g ^f(p, c), 
c(x, t) = cb > as \x\ — > oo 
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where c denotes the density of nutrients, and cb the far field supply of nutrients (from blood 
vessels). The coupling functions $, $ are assumed to be smooth and to satisfy the natural 
hypotheses 

<9 P $<0, <9 C $>0, $( PM)Cb ) = 0, 

(1.5) 

d p ty<0, d c ^>0, *(p,0)=0. 

Variants are possible; for instance, we could assume that nutrients are released continuously 
from a vasculature, then leading to an equation as 

d t c — Ac = —q *$?(p, c) + r(ce — c). 



We will consider below the purely mechanical model (1.1) under assumptions ( 1.2 )— ( 1.3 ), 



and also the system (1.4), where nutrients are taken into account, under assumptions (1.3) 



and (1.5). In both cases we will show that the asymptotic limit m — > oo yields a free 



boundary model of Hele-Shaw type, as we were looking for. 



Let us recall that the mathematical theory of the limit for an equation like (1.1) in the 
absence of a growth term is well developed, the asymptotic limit being in this case the 
standard Hele-Shaw model for incompressible fluids with free boundaries. Early papers on 
the subject 0, El IH1 El 1221 [2H [23, EH] consider situations in which mass is conserved and 
the limit is stationary. In order to have a non-trivial limit evolution one needs some source, 
either in the equation or at the boundary of the domain. A first example, in which there 
is an inwards flux at infinity, is given in [2], where the authors study the limit m — > oo 
of self-similar focusing (hole-filling) solutions to the porous medium equation u t = Au m . A 
second example, in which there is a bounded boundary with a nontrivial boundary data, was 
first considered in [28], and later on in [221 E2J EH ES]- To our knowledge, the present paper 
is the first one in which the evolution in the limit is produced by a nonlinear source term in 
the equation. In order to pass to the limit, three approaches have been used: weak solutions, 
variational formulations (using the so-called Baiocchi variable), and viscosity solutions; see 
[3"3"| El] for this last case. The weak formulation of Hele-Shaw was first introduced in [20J, 
and the variational formulation in [23J. 

Let us also mention that the Hele-Shaw graph, and hence the Hele-Shaw equation, can be 
approximated in other ways, for example by the Stefan problem. This situation has also 
been considered in the literature, even in situations where there is an evolution in the limit; 
see for instance [TQ1 EH El] ■ 

We would like to stress that including the growth term is not a simple change: several 
powerful but specific tools do not apply any longer. More deeply, as we explain later, several 
approaches have no chance to work as they are. This is the reason why we will work on 
an equation for the cell density itself rather than the pressure. We also would like to point 
out that models of Hele-Shaw type are still an active field arising in several unexpected 
applications, see for instance [2T], and that surface tension is not covered by the present 
work; see [igilEl 021 EE]. 

Organization of the paper. We first study the simplest mechanical model in Section [2} The 
more complex model with nutrients is studied in Section |4} The uniqueness proofs for both 
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cases are highly technical, and are performed separately, in sections [3] and [5] These two sec- 
tions can safely be skipped by the readers who are mainly interested in the applications. We 
finally include an appendix devoted to some interesting examples for the purely mechanical 
model, that illustrate several phenomena. 

Some notations. We will use several times the abridged notations g(t), pit), meaning 
g(t)(x) = g(x,t), p(t)(x) = p(x,t). Given any T > 0, we denote Qt = R N x (0,T), while 
Q = R N x (0,oo). 



2 Purely fluid mechanical model 



We start our analysis by a detailed investigation of the purely mechanical model (1.1 ), which 



does not take into account the consumption of nutrients, with a pressure field given by (1.3). 



A simple change of scale allows us to assume without loss of generality that g c — 1. We 
arrive to the porous medium type equation, set on Q, 



(2.1) 



d t g = Ag m + g$(p), p = P m {g) 



m 



jm—X 



-Q 



m — 1 

We assume that the initial data g^ are such that, for some g° E L+(M. N 

Qrn > °) P m(.Qm) < PM, 

0, \\^x ^ g ^ L\\I^mN^ < C, i 



g(0) = g° m . 



(2.2) 



\Qm ~ Q 



1,...,N. 



2.1 Main results 



Free boundary limit. Let (g m ,p m ) be the unique bounded weak solution to (2.1). We 



will prove that, along some subsequence, there is a limit as m — > oo which turns out to be a 
solution to a free boundary problem of Hele-Shaw type. 



Theorem 2.1 Let $ and {f?„} satisfy (1.2) and (2.2) respectively. Then, after extraction 
of subsequences, both the density g m and the pressure p m converge for all T > strongly in 
^(Qt) asm-)- oo to limits g^ E C([0, oo); L 1 (R N )) nBV(Q T ), p^ E BV{Q T ), that satisfy 
< £oo < 1, < Poo < p M , and 



(2.3) 



d t g 00 = A Poo + g 00 $(p tx> ) inV'{Q), eoo (0) = g° in L 1 (R N ) 



plus the relation p^ E -Poo(^oo); where is the Hele-Shaw monotone graph 



(2.4) 



Poo(g) 



o, 

[0,oo) 



< Q < 1, 
= 1. 



Note that (2.4) means that a.e. P^ > and P^ = a.e. in {0 < g < 1}. 
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We will obtain also several important qualitative properties for the limit. On the one hand, 
the 'tumor is growing' and the pressure increases, that is, 

(2.5) d t6oo >0, d tPoo >0 mV'{Q). 

On the other hand, if the initial data {p^} have a common compact support, then the limit 
solution propagates with a finite speed: Qooit) and Poo(t) are compactly supported for all 
t > 0. Finally, the constructed limit solution enjoys the monotonicity property, inherited 
from the case where m is finite, 



Qoo{t) > Qoo(t), 
Poo(t) > Poo(t). 



Transport equation and complementarity formula. We will also obtain L 2 esti- 
mates on the gradients Vp m that will show on the one hand that 

Vpoo G L 2 (Q T ) for all T > 0, 

and on the other hand that (Qoo,Poo) solves a transport equation. 



Theorem 2.2 The limit solution (Qoo,Poo) obtained in Theorem 2.1 satisfies 

d t Qoo ~ div (^ooVpoo) = Qoc®(Poo) 

in a weak sense. 

A direct calculation shows that the pressure p m satisfies 

(2.6) dtp m = {m- l)pmAp m + \Vp m \ 2 + (m - l)p m ${p m ). 
Hence, if we let m — > oo, we formally obtain the complementarity formula 

(2.7) p 00 (Ap 00 + $(p 00 )) =0. 
This formal computation can be made rigorous. 

Theorem 2.3 Under the assumptions of Theorem 2.1, the limit pressure p^ satisfies the 
complementarity formula 

(2.8) / (— |Vpoo| 2 + Poc^iPocj) = for almost every t > 0. 
Jr n 



It is worth noticing that (2.8) is equivalent to the strong convergence of Vp m in L (Qt) for 
all T > 0; see Lemma 2.5 



Since the limit pressure p^ is expected to be continuous in space for all positive times, the 
positivity set 

fi(t) := {x; Poo(x,t) > 0}, 
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should be well denned for all t > 0. Notice that it coincides almost everywhere with the set 
where g^ — 1; see Figure [Tj Indeed, on the one hand, by the definition of the graph we 
have Q(t) C {x; Qoc(x,t) = 1}; on the other hand, if we had p^ = and g^ — 1 in some 
set with positive measure, then would continue to grow (exponentially) there, which is 
a contradiction. Therefore, Q(t) may be regarded as the tumor, while the regions where 
< Qoo < 1 (mushy regions, in the literature of phase-changes) correspond to precancer 
cells. 





Figure 1: Effect ofm large. A solution to the mechanical model (1.1 ), (1.3) in one dimension 
with $(p) = 5(1 — p). Left: m = 5. Right: m = 40. The upper line is £>; the bottom line is 
p (scale enlarged for visibility). Notice that the density scales are not the same in the two 
figures. The initial data is taken with compact support and the solution is displayed for a 
time large enough (see Figure [2] below for and intermediate regime). 



The complementarity formula (2.7) indicates that the limit pressure at time t, Poo(t), should 
solve the elliptic equation 



(2.9) 



A Poo {t) = $(poo(*)) in tl(t), poo(t) e H%(n{t)). 



a problem which is wellposed if Q(t) is smooth enough, since pi—)- $(p) is decreasing. This 
implies in particular that in general regularity in time for the pressure is missing, since time 



discontinuities may show up when two tumors meet; see Subsection A. 2 in the Appendix. 



Geometric motion vs. equation on the cell number density. To complete the 
description of the limit problem we should be able to trace Q(t) starting from its initial 



position. The pressure equation (2.6) suggests that we should have 

dtPoc = |Vpoo| 2 at dQ(jt), 
which leads to a geometric motion with normal velocity V at the boundary of Q(t) given by 
(2.10) V=\V Poo \. 

Thus we have arrived to a geometric Hele-Shaw type problem, which is the classical one 
when $ = 0. 

The above formal computation is expected to be true if we prescribe a fixed initial pressure 
Pm(0) = p° (which implies that the initial densities converge to the indicator function of 
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the positivity set of p°). When $ = 0, it was proved to be true in a viscosity sense in 
see also [28J for an earlier result in this direction using a variational formulation of Baiocchi 
type. 

However, if the initial densities g^ are such that g° is below 1 in a set with positive measure, 
the result is no longer true. The main point is that the tumor may meet precancer zones. 



At a meeting point (x,t), the example in Subsection A. 3 in the Appendix suggests that the 
tumor grows faster (also for the case $ = 0), with a normal velocity given by a rule of the 
form 



(2.11) 



V 



|VPoo| 



where g is some limit of g as t — > t and x — > x from the outside of Q(t). We leave open this 
problem, which seems to be a challenging extension of the viscosity method in 



Even if we consider the modified geometric motion law (2.11 ) instead of (2.10), the geomet 



ric formulation does not carry all the information of the limit solution. Indeed, the density 
in precancer zones evolves, with an exponential growth, until it reaches the level g^ = 1, a 
fact that is not captured neither by ( |2.9|), n or by (2.11); see Figure [2] and the examples in the 
Appendix. When $ = we also need (12.3 ) to give a full description of the limit. However, 



in this latter case the evolution in the mushy regions is simpler, since the density does not 
change until it interacts with the positivity set for the pressure. 



Figure 2: Cell density and pressure carry different informations. Here m = 40 and the 
initial data g is less than 1. The solution is displayed at four different times. It shows how 
the smooth part of g strictly less than 1 is growing with p = (figure on the left). When g 
reaches the value 1, the pressure becomes positive, increases and creates a moving front that 
delimitates the growing domain where g ~ 1. Thin line is g and thick line is p as functions 
of x. See also Figure [3] in the Appendix. 

Let us point out that, in the case of a compactly supported initial data, precancer zones 
disappear in finite time. 

Uniqueness. To complete the description of the asymptotic limit, we conclude with a 



uniqueness result, relying on a Hilbert duality method, for the free boundary problem (2.3) 



Theorem 2.4 There is a unique pair (g,p) of functions in L°°((0, oo); L 1 (IR Af ) fl L c 
g G C([0, oo); L 1 (IR Ar )) , p G Poo(q), satisfying (2.3) and such that for all T > 0: 



p{t) is uniformly compactly supported for t G [0, T]; 
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(ii) I Vp| G L 2 {Qt) 

Ciii) d t peM l (Q T ), dtgeM^Qr). 



It is important to notice that if, in addition to ( |2.2[ ), we assume that the initial data 

9° 
- m 



^ are uniformly compactly supported, then the limit solution (^oo,Poo) to (2.3) given by 



Theorem 2. lj falls within the uniqueness class. As a consequence, convergence is not restricted 
to a subsequence. 



Remark. The compact support assumption in Theorem |2.4 can be removed, to the cost of 
using rather technical arguments, in the style of the ones employed in the uniqueness proof 



for the case with nutrients, Theorem 4.2 below. We keep it for the sake of clarity. 



The rest of this section is devoted to the proofs of the above mentioned results, except 
uniqueness, which is postponed to Section [3j 



2.2 Estimates, strong limit and existence 



We consider the solutions to ( 2.1[ ) and proceed to obtain the estimates that allow us to pass 
to the limit and prove Theorem 2A_ This is a standard porous media type equation and 
all the manipulations below can be justified, see the monograph [42J. The main observation 
here is that the lower estimate on Ap m that does the job when passing to the limit in porous 
media equations without a source, does not work in the present case. However, a control 
can be obtained on the quantities Ap m + $(p m ), which will be enough for our purposes. 

L°° bounds FOR Q m , p m . Standard comparison arguments yield 



< Qm < 



m 



1 



m 



-Pm 



l/(m-l) 



< p m = P m {Qm) < Pm- 



This allows in particular to avoid initial layers (which are present when 



> 1; see 



)■ 



L 1 bounds FOR Q m , p m . Let g m , Qm be two (non-negative) solutions to (2.1), and 
Pm = Pm(Qm), Pm = Pm{.Qm) the corresponding pressures. We have 

d 



dt 



{Qmit) - Qmit)}- 



< 



+ 



®(Pm(t)){Qm(t) - Qm(t)} + 

Q m (t)($(p m (t)) - $(]3 m (t))) sign + (£ m (t) - Q m (t)) 
< $(0) / {Qm{t) ~ Qm{t)}+, 



{Qm{t) ~ Q m {t)}+ < e* (0) * / {Qm(0) - Q m (0)} + - 



from where we get 



This "almost contraction" property yields the uniform (in m) bound 



(2.12) 



<e mt hl\W 



< Ce^ 0)t . 
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On the other hand, p n 



m—1 1 



' m—1 

> m 



Pr, 



!. Thus, using (2.12) we conclude 



lbm(*)IUi( R ")<Ce*^ form>2. 



A semiconvexity estimate FOR p m . According to the hypotheses (1.2) on the growth 
function, r$ = min ($(p) — pQ'{p)) > 0. We will prove that 

P&[0,Pm] 

(2.13) Ap m (t) + ®(p m (t)) > -r e-( m - 1 >*V(l - e"^" 1 ^*'). 

As a consequence, in the limit we will have Apoo + $(poo) > 0. 



We borrow the idea to prove inequality (2.13) from [TJ, where the case $ = was considered. 



"m) j 



We write the equation (2.6) as 

(2.14) d t pm = {m - l)p m w + |Vp m | 2 , with w = Ap m + $(p m ). 

Let us denote t> = Ap m . Since $' < 0, we have 

d t v = (m - l)p m A«j + 2(m - l)Vp m ■ Vw + (m - l)vw + 2Vp m -Vv + 2 Y^iji^XiXjPr, 

dt($(pm)) = &(p m )d t Pm = (m- l)p m $'(p m )w + &(p m )\Vp m \ 2 
> (m - l)p m $'(p m )w + 2V($(p m )) ■ Vp m , 

which gives 

<9 t u> > (m - l)p m Aw + 2mVp m ■ Vw + (m - l)w 2 - (m - l)($(p m ) - p m $'(p m ))w. 
The function = — r,j>e _ ( m ~ 1 ) r **/(l — g-C" 1 - 1 )*"**) is a subsolution to this equation, 



and (2.13) follows. 



Remarks, (i) The right hand side of (2.13) behaves as —\/{t(m — 1)) for t ~ 0. 

(ii) For $ k we have r$ 0, and we recover the well-known result for the case $ = 0, 
namely Ap m (t) > —l/(t(m — 1)). 



Bounds for <9 t p m , <9t£ m . We combine (2.13) with (2.14) to obtain an estimate from below 
for the time derivative of the pressure, 



(2.15) 



d t Pm(t) > -(m - l)p m {t)r^- 



for t > 0. 



J g— (m— l)r$i 

This in turn gives an estimate from below for the time derivative of the density, 

g— (to— l)r $ t 



(2.16) 



g— (m— l)r$t ' 



m — )■ oo in (2.16) and (2.15). 



The monotonicity inequalities (2.5) for the limit problem are then obtained just by letting 
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We now use \d t Q m \ = d t Qm + 2{<9 i £ m }_ to obtain 
\d t g m (t)\\m^) = ~n / Qm(t)+2 / {d tQm {t)}_ < $(0) + " " 



dt L N " mv y ' /icpjv ^™^ J - - V v 7 ' 1 - e-("-iM 



This, together with (2.12), leads to a uniform bound in space and time for time intervals of 
the form te [^,T[. 



An analogous computation shows that 

/ / \dtPm\ < / Pm(T) - / 

Pm I -. I 

JiJrn J r n J r n \m-lj 



m— 1 



T / e —(m— l)r*t 



m — 1 



+(m-l)r« / ( t _ (n ^ 1)r , t / ) dt < C(T). 



L 1 bounds FOR V^ m , Vp m - Let a = min pe [ 0jPM ] |$'(p)| > 0. We consider the equation 
for d Xi g m , multiply it by siga(d Xi g m ) = sign(<9 x .p m ) , and use Kato's inequality; thanks to the 
monotonicity of $ we obtain 

d t \d x .Q m \ - A(mg™~ 1 \d Xi g m \) < $(p m )\d Xi Qm\ + \Pm) Q m \d Xi p m \ 

iQm\ QQm 

Integrating in Q t , we get 

\\d Xi e m (t)\\m^) + otjj Q m \d XiPm \ < \\d Xi g m \W { R»)e mt < Ke*®*, 



which yields, on the one hand, that 

and on the other hand that 

\\d Xi p m \W{Q T ) < // rng™- 2 \d Xi Q m \+2 / / Q m \d Xi p m \ < C(T). 

J jQ T n{ Sm <l/2} J JQ T n{g m >l/2} 

Convergence and identification of the limit. Since the families g m and p m are 
bounded in W lo ' c (Q), we have strong convergence in L\ oc (Q) both for g m and p m . 

To pass from local convergence to global convergence in L 1 (Qt), we need to prove that the 
mass in an initial strip t e [0, 1/R] and in the tails \x\ > R are uniformly (in m) small if 
R is large enough. The control on the initial strip is immediate using our uniform, in m 
and t, bounds for ||^ m (t)||z,i(MJV) and UPmCOIIi 1 ^)- ^ n OT( ^ eY to control the tails, we consider 
ip e C°°(R N ) such that < tp < 1, (p(x) = for \x\ < R - 1 and (p(x) = 1 for \x\ > R, and 
define <£r(x) = cp(x/R). Then, for any m > 2, 

4 / Qm(t)ip R < 2R' 2 \\g m (t)\\^} RN) \\g m (t)\\ L i {R N ) \\Aip\\ L ^^N ) + ^(0) [ Q m (t)<pi 
dt Jrn *• > J R N 

< CR- 2 + $(0) / Q m (t)<p R . 
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Therefore, for any t G [0, T] we have 

0< [ g m (t)<p R < e*w([ Q m y R + CR- 
Jrn \Jrn 

< e* (0)T ( \\qI - Q°\\ L i iRN) + [ g° VR + CR- 2 t) < e 
V Jr n J 

for R and m are large enough. Since 

/ Pm(t)p R <2\\Q m (t)\\™-} N) / g^fl 

for any m > 2, the tail control for the pressures follows easily. 

After extraction of subsequences, we can pass to the a.e.-limit in the equation g m Pm = 
Pm +m ^ m to obtain that (1 — Qoo)Poo — 0. Passing to the limit in the estimates, we have also 
< £oo < 1, < Poo < p M and Qoo,Poo G BV(Q T ) for all T > 0. 



All the above is enough to prove that the pair (f?oo)Poo) solves the Hele-Shaw equation (2.3) 
but for the question of the initial data. 

Time continuity and initial trace. We need delicate arguments that are fortunately 
in the folklore of the topic of nonlinear diffusion equations, and borrow also the conclusion 
from |3Qj. Because is non-decreasing in time, we can write for a test function < ((x) < 1 
and < ti < t 2 < T, 

/ |0oo(*2)-0oo(tl)|C = / [goofa) - 0oo(*O] C = / / [P00AC + £oo$G0oo)C] 

Jr n Jr n Jti Jr n 

< C(T) (tj, - tO ( || Aflloo + 1). 

Taking a sequence of such uniformly smooth functions that converge to 1, we find that 
goo G C([0, oo); L 1 ^^)) (in fact with a locally uniform Lipschitz constant). 

In order to identify the initial trace, we observe that for any test function ( as above, 

/ Qm{t)(- [ Q°mC= [ [ (p m AC + £? m $(p m ))C 

Jr n Jr n Jo Jr n 

Letting m — > oo, we have 

/ Qoo(t)(- I g° 00 (= [ [ (PooAC + £?oo$(Poo))C 

Jr n Jr n Jo Jr n 

Letting first t — > and then ( — > 1, we conclude that g(0) = g° in L 1 (IR Ar ). 

2.3 L 2 estimates for Vp m , transport equation and complementarity 
formula 



We now prove Theorems |2.2| and |2.3[ both of them connected with the fact that Vp m G 
L 2 (Q T ) for all T > 0. 
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Since we have already proved that both g m and p m converge strongly in L 1 (Qt) for all T > 0, 
the first of these theorems just depends on obtaining a uniform bound in L 2 (Qt) for Vp m . 



Proof of Theorem 2.2. We rewrite the pressure equation (2.6) as 

1 



(2.17) 



dtPr, 



m 



Ap 2 m - (m - 2)\Vp m \ 2 + (m- l)p m $(p r . 



2 



Integrating in Qt we obtain the required estimate, 



|| Vp 



m\\L 2 (Q T ) 



< 



^(0)|bm||Li(Q T ) + 



m 



m II L 1 



□ 



Remark. Combining (|2.17|) with (|2.15|), we obtain 

1 



l|Vp„ 



< 



Pm(t) 



III 



m — 2 



$(0) + (to - l)r$ 



J g— (m— l)r$t 



for all m > 2. This gives a uniform estimate for ||Vp m (t)||£,2( K jv), t G [t,T], for any fixed 
values < t < T. However, it does not allow to go down to t — 0. 



From the equation in (2.1) and the L l bound on dtQ m we conclude that 

||A^(t)|| Llf ^<Ce*(°)*. 



This is an optimal bound since p^ has 'corners' on dQ(t). Because this only gives space 



compactness, it is not enough to establish the complementarity formula (2.7) on the pressure 



when passing to the limit in (2.6). We will need to perform a time regularization argument 
a la Steklov. 



Proof of Theorem 2.3, Let u e (t) > be a regularizing kernel in time with support (— e, 0). 



We use the notation g m)E (t) = u £ * g m (t) = J R w e (t — s)g m (s) ds. The equation (2.1) gives 
(2.18) A(g™*u e ) = d t g mj£ - (g m $(p m )) * u e . 

From here it follows that U m = A(^™ * u> £ ) is uniformly (in m) smooth in time and H 1 
in space for e fixed. Therefore U m converges strongly in L\ oc (Q). Hence, after multiplying 



(2.18) by p m , we may pass to the limit to obtain 

PaofAfyoo * U e ) + (^oo^(Poo)) *^ £ ) = lim Pmd t 0m,e- 

m— >oo 

In order to estimate the right hand side we make the following decomposition, 



(p m dtQm, £ )(t) 



rn 



m 



QZ (s)d t g m (s)ui £ (t - s)ds 



lm(t) 



+ 



ni 



m — 1 



C 



(eZ-\t) - gZ~\s)) d t g m (s) + - u £ {t - s)d 



Cm 



limit) 



m 



v 

IHm(i) 
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where C is a constant such that d t Q m + 7 > 0; see estimate (2.16). 
For the first term we have 

1 f +E / . \ff<^ m ,\<\, . C, 



lm(t)\ < 7 / -«) / l3t*C00l < 



m — 1 L \ mN 1 m — 1 



Regarding the second term, the estimate (2.15) implies that d t p m > —C. Therefore, 
£m _1 (^) — £m _1 ( s ) — Ce. Let C G V(Q), and r the smallest time in its support. Then, 

C Hm < J J J C ^£m( s ) + —^j u e {t — s) dsdxdt < Ce. 
The third term is easy to treat. Since s > t > 0, for any test function ( as above, 

// CHI W — ► -C // C / (Poc(t) - Poc(a)) ^ ~ g) ttotfadt = 
JJq m -+°° JJq Jr s 

as e — > 0. 

These three calculations give, in the limit m — > 00, 

Poo (A^ ★ o; e ) + (g 00 ^(p 00 )) *u E ) < o(l) as £ ->■ in 

To recover the desired information, 

p 00 (Ap 00 + $(poo)) < 0, 

it remains to pass to the limit as e — > after noticing that PooQoo = Poo and that the 
differential term can be treated through its weak formulation; indeed, after testing it against 
a test function as before, it is written as 



Q 



[CVpoo • VOoo * U £ ) + P00VC • V(poo * U) £ )] 

which passes to the limit e — > because we already know that Vpoo G L 2 (Qr)- Therefore, 

~ Jj [C|Vpoo| 2 +PooVC- Vpoo + CPoo$(Poo)] =0. 



□ 



As we have already mentioned, the complementarity formula is equivalent to the strong 
convergence of the gradients. 



Lemma 2.5 The complementarity formula (2.8) holds if and only if Vp m — > Vpoo strongly 

' m— >oo 

m L 2 (Qt) /or a// T > 0, 



13 



Proof. We consider a test function ( G T)'{Q) and use (2.14) and (2.3) to obtain 



(p m ( A Pm + $(Pm)) = — T IIC (dtPm ~ |Vp m | 2 ) > 0. 



Q m ~ 1 



This implies, after an integration by parts, that 

// [-C|Vp m | 2 -p m VC- Vp m +Pm(®(Pm)] > 0. 

J JQ m ^°° 

Let r = liminf m _ 5 . 00 ffqCl^Pml 2 — JJq CI Vpoo| 2 - Thanks to Fatou's Lemma, r > 0. On the 
other hand, since p m converges strongly, after extraction of subsequences we have 



and r = is equivalent to the strong convergence. □ 
2.4 Finite speed of propagation 

If the initial data are compactly supported uniformly in m, we can control the supports of 
Pmit) and Q m (t) uniformly in m and t G [0, T] for all T > 0. This implies that the speed of 
propagation is finite for the limit problem. 

The control is performed comparing the pressures with functions of the form 



(2.19) p (l , t)= ^__M 2 () 

The key point is that these functions, which are viscosity solutions of the Hamilton- Jacobi 
equation P t = |VP| 2 , are supersolutions to the equation 

(2.20) d t p = (m - l)pAp + | Vp\ 2 + (m - l)p$(0) 

for some time interval which does not depend on m. This idea was already used for the case 
$ = in [2EJ I2H]- However, in the absence of reaction, P is a supersolution for all times, 
and the proof is a bit easier. 



Lemma 2.6 Let {g^} satisfying (2.2) and such that all their supports are contained in 
a common ball Br- Then, for every T > 0, there is a radius TZt depending only on 
sup m ||^||l°°( K jv), R and T, such that the supports of {g m {t)} are contained in the ball 
B-Rt f or all m > 1 and t G [0, T] . 



Proof. We consider P as in (2.19) with r = A^/(4$(0)) and C large enough so that 



P(x,0) > 

Qm( x ) f° r ^1 x ^ 1R and m > 1. An easy computation shows that 

' N 
Xr + t) 



d t P - (m - 1)PAP - | VP| 2 - (m - 1)P$(0) = (m - l)P(x, t) ( — - $(0) I > 
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for all t G [0, 



N 



||1(0) . Since the functions p m are subsolutions to (2.20), we conclude that the 
supports of the functions g m are contained in the ball of radius ^/i^jjy f° r ah t G [0, ^p-]- 
Moreover, UPmC ^irnO lU 00 ^) — ^> an< ^ the argument can be iterated to reach the time T in 



^4*(0). 

a finite number of steps. 



□ 



3 Uniqueness for the purely mechanical limit problem 



In this section we prove, under suitable assumptions, that the limit problem (2.3) has a 



unique solution, Theorem |2.4[ The main difficulty comes from the fact that p is not a 
Lipschitz, single-valued function of g. Hence, we cannot apply directly the ideas developed 
in [7] to adapt Hilbert's duality method to the porous medium equation. The technique has 
to be modified, in the spirit of |17j . 

Dual problem. Consider two solutions, (qi,Pi), (£2,^2)- Let fl be a bounded domain 
containing the supports of both solutions for all t G [0, T], and Qt = x (0, T). Then we 
have 

(3.1) ff [(Qi-Q2)d^ + (p 1 -p 2 )A^ + (Q 1 $(p 1 )-Q 2 $(p 2 ))^] =0 
for all suitable test functions ip. This can be rewritten as 

(3.2) ff (g 1 -g 2 +p 1 -p 2 )[Ad t ^ + BA^ + A^(p 1 )^-CB^] = 0, 



J Jn T 

where, for some fixed v > 0, 

0<A = - f 

[Ql ~ Q2) + [Pi - P2> 

0<B = - ( ,<1, 

{Qi - 02) + [Pi ~ P2> 

o< c = <„. 

Pi -P2 

To arrive to these bounds on A and B, we define A = when Qi = g 2 , even when p\ = p 2 , 
and B = when pi = P2, even when Q\ = q 2 . 

The idea of Hilbert's duality method is to solve the dual problem 

Ad t ip + BAip + A<$>( Pl )ip - CBi) = AG in Q T , 
ip = in dfl x (0, T), ip(-,T) = in fi, 
for any smooth function G, and use ^ as test function. This would yield 

(Qi ~ Q2)G = 0, 



T 
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from where uniqueness for the density is immediate. Uniqueness for the pressure would then 



follow from (3.1 ). 



However, on the one hand the coefficients of the dual problem are not smooth, and, on the 
other, since A and B are not strictly positive, the dual equation is not uniformly parabolic. 
A smoothing argument is required. 

Regularized dual problem. Let {A n }, {B n }, {C n }, {$i, n } be sequences of smooth 
bounded functions such that 



\\A - A n \\ L 2 {nT) < a/n, 
\\B - B n \\ L 2 {QT) < 0/n, 
\\C - C n || L 2(Q T) < 5i/n, 
k ||$(pi) - $i,n|U2(n T ) < S 2 /n, 



l/n < A n < 1, 
1/n <B n <l, 
0<C n < Afi, 

|$l,n| < M 2 , 



\dtC n \\ L i^n T ) < Ki 
|V$i,„|| L 2 ( n T) < K 2 , 



for some constants a, (3, 61,62, Mi, M 2 , K\,K 2 . 
For any smooth function G, we consider the solution ip n of the backward heat equation 



(3.3) 



B B 

d t ipn + -j^A^n + $l, n ^n - C n -^^ n = G in Q T , 

ij n = mdttx(0,T), ^ n (-,T) = in a 



The coefficient B n /A n is continuous, positive and bounded below away from zero. Thus, the 
equation satisfied by ip n is uniformly parabolic in Qj<. Hence ifi n is smooth and can be used 



as a test function in (3.2) 



Combining (3.2) and (3.3), we have 

// { ei -Q 2 )G = 11-11-11 + 1, 

where 



4 

ni 





In T 


((ft 


- 02) + (pi 


-p 2 ))^(A-A n ) (A^ n -C7 n ^n 




In T 


((ft 


- 02) + (pi 


-p 2 )){B-B n ) (Aif) n -C n ip n ), 




In T 


(ft ~ 


ftO($(pi) " 






In T 


((ft 


- fe) + (pi 


-p 2 ))B(C-C n )i> n . 



Limit n. ->■ 00. Uniform bounds and conclusion. Our aim is to prove that 
lim^oo /* = 0, i = 1, ... ,4, which implies that Qi = g 2 . For this task we will use some 
bounds that we gather in the following lemma. 
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Lemma 3.1 There are constants K{, i = 1,2,3, depending on T and G, but not on n, such 
that 

\\i>n\\L°°(ClT) ^ ^1) SU P ||Wn(*)|U a (n) < K 2, \\{B n /A n ) l/2 (Alp n - Cntpn) ||l2(H t ) < «3- 
0<t<T 

Proof. The first bound is just the maximum principle because C n is non-negative and 
uniformly bounded. 



The second bound is obtained multiplying the equation in (3.3) by Aip n — C n ip n - After 
integration in Q x (t,T), we get 
(3.4) 



B 



t Jn Ai 



-\Alp n - C n tp n \ 



(t) 



+ 

>t Jn 
< K ( 1 — t 



T 

l|V^(s)|||2 (n) ds 

where we denote by K various constants independent of n. We now use the Gronwall lemma 

to 
□ 



to conclude our second uniform bound. Then we can use again the inequality (3.4) to 
conclude the third one. 



We now get 



\ll\ < Kll -^\A-A n \\A^ n -CM </q(£ n /^) 1/2 (A-A n )|U 2(fV) 



,1/2. 



< Kn^WA- A n \\ L 2 (nT) < Ka/n 

\I 2 n \ < K [[ \B-B n \ \A^ n - CM < K\\{A n n l l 2 lB n fl\B - B n )\\ L , {n 
J Jn T 

< Kn 1/2 \\B - B n \\ 2 < Kf3/n 1/2 ; 
\qi - g 2 \ |$(pi) - $i, n | \ipn\ < K\\$(pi) - $i,n|U2 (nT) < K5 2 /n; 



t) 



\r\ < k 



n T 



B\C - C n \ |^| < K\\C - C n \\ L 2 [nT) < K/n. 



Once we have proved that Qi = Qz, equation (3.1) says that 



// [(pi-P2)A^ + gi($(pi)-$(p s 
J Jn T 



from where the result will follow by taking as test function if) = pi—p 2 , using the monotonicity 
of $. 
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4 Model with nutrients 



When nutrients are taken into account, we are led to problem (1.4). We again assume that 



the pressure field is given by (1.3), with g c set equal to 1. 



Our assumptions on the initial data are stronger than for the purely mechanical problem. 



Namely, in addition to (2.2), we assume that for some c° such that cb — c° G L+(R ), 



< C° < C B , 



c Wl 1 



hiWL 1 



<c, 



1,...,N, 



pv(<&Vp£J + Q°m Hp° m ,c J\\ LHRN) < C, \\Ac° m - el *(p° m ,c m )\\ LHRN] < C 



(4.1) 



An interesting difference with the purely mechanical model is that when nutrients are not 
enough the cells might die (or become quiescent), which is represented by $(p, c) < for c 
small enough with p given. These cells, which are in principle in the center of the tumor, 
are replaced by cells from the boundary moving inwards by pressure forces. In particular 



the growth inequalities (2.5) cannot hold true here. 



A typical choice is ^(c,p) = c, $(c, p) = $(p)(c + Ci) — c 2 with q > two constants and 
$(p) a function having the same properties as in the purely fluid mechanical model. 



4.1 Main results 

As in the case without nutrients, in the limit we obtain a free boundary problem of Hele-Shaw 
type. 



Theorem 4.1 Let®, ^ satisfy (1.5), and {(g m ,c Q m )} satisfy the hypotheses (2.2) and (4.1) 



Then, after extraction of subsequences, the density g m , the nutrient c m and the pressure p m 
converge for all T > strongly in L 1 {Qt) as m —¥ oo to limits f?oo,Coo,Poo G BV(Qt) that 
satisfy < g^ < 1, < < cb and < poo < pu, and 



(4.2) 



dtQoc = Ap^ + ^oo$(Poo, Coo), £oo(0) = Q° , 
dtCoo = ACoo - ^oo^(Poo, Coo) Coo(0) = C°, 

in a distributional sense, plus the relation p^ G -Poo(f?oo); where is the Hele-Shaw 



graph (2.4). 



When we say that (^oo,Poo, Coo) is a solution to (4.2) in a distributional sense, we mean that 



// {goo4>t +PooA^ + Qoo^iPocCoo)^} = ~ / g°lj)(-,0), 

JJq Jrn 

{CooCt + CooAC - 0oo*(Poo, Coo)C} = - / C°C(-, 0), 

Jrn 



for all test functions tpX G C£°(Q). 
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As mentioned above, for the system we lose the monotonicity properties (2.5), and we are 
not able to prove the expected continuity of g^ and cb — Cqo in L 1 (IR Af ). However, we are 
still able to prove that Vpoc in L 2 {Q T ). Therefore, the equation on the cell density can also 
be written transport equation, 

(4.3) <9t£oo - div (^ooVpoo) = £oo$(Poo, Coo). 

If the initial data g^ are compactly supported uniformly in m, we can control the supports 
of g m (t) uniformly in m for each t > 0. Hence, in the limit problem tumors propagate with 
a finite speed (this is not true for the nutrients), and a free boundary shows up. This can be 
proved through a comparison argument with the same barrier functions as in the scalar case. 
Moreover, the constructed solution will fall within a class for which we prove uniqueness. 

Theorem 4.2 There is a unique triple (g,p, c), g,p, cb — c G L°°((0, oo); L 1 (IR iV ) nL°°(]R JV )) 7 



P £ Poo(q), satisfying (4.2) in the distributional sense and such that for all T > 0: 

(i) p(t) is uniformly compactly supported for t G [0, T]; 
(it) |Vc|,|Vp| eL 2 (Q T ); 
(in) d t p G M\Q T ), d t g G M\Q T ). 

As in the purely mechanical model, we can write an equation for the pressure p m , 

(4.4) d t p m = (m - l)p m Ap m + I Vp m | 2 + (m- l)p m $(p m , c m ). 
This suggests that the complementarity formula 

(4.5) Poo [A Poo + $(p )) = o 

holds. However, we have not been able to establish it rigorously. The reason is that, in 
contrast with the case without nutrients, we have failed to control Ap m + $(p m ,c m ) from 
below by means of a comparison argument, or to prove the strong convergence of Vp m in 
L 2 {Qt) for all T > 0. We leave open the question to give conditions on $(p, c) allowing to 



prove formula (4.5). 



The rest of this section is devoted to prove the first of these theorems. The uniqueness 
result is postponed to a later section. 

4.2 Estimates 



L°° BOUNDS FOR £> m , p m , c m . The assumptions on the growth functions (1.5) imply, using 
standard comparison arguments, 

< g m < [ Pm < p m {x, t) < p M , < c m < c B . 

\ m I 
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L 1 bounds FOR g m , p m , c m . Let (g, c) and (g, c) be two solutions of the system for a fixed 
m, with corresponding pressures p and p. We subtract the equation for g from the equation 
for g, and multiply the resulting equation by sign(f) — g), and do an analogous manipulation 
for the c variable. We integrate by parts, and obtain 

(e®(p, c) - Q$(p, c)) sign(j» - p) (t), 
(q&(p, c) - g^(p, c)) sign(c - c) (t). 




\\(q~ Q)(t)\\ L i(R») < 

|| (C-C)(t) || L l (R iY) < - 



Let jj, be a positive constant to be chosen later. We have 



j t (|| (e - Q)(t) IUhrat) + /i|| (c - c)(f) lUi(RiV)) < 

/ (q - q) ($(p, c) s\gn(g - g) - fi9(p, c) sign(c - c)) (t) 
Jr^ l j 

V v ' 

X 

+ / £(($(p, c) - $(p, c)) sign(e - g) - ji(*(p, c) - *(p, c)) sign(c - c)) (i) 

JRN L J 

N v ' 

+ / q{(®(p, c) - $(p, c)) sign(^ - £) - n(V(p, c) - *(p, c)) sign(c - c)) (f) 



K 



Let A = {0 < p < p M , < c < c B }, a = mm A \d p $\ > 0, (3 = m&x A \d p ^\. Since <9 P $ < 0, if 
< fj, < a/ (3, there are constants C independent of m such that 



X< 



\g- g\($(p,c) + nV{p,c)) (t) < C\\(g - g)(t)\\ L i {RN) , 



J< \\p\\l~(Q) {\\dM L ~ { A)+PL\\dML°°(A)) \\(c-c)(t)\\ LHRN) <C\\(c-c)(t)\\ Ll 

K < f \g{- |$(p,c) - $(p,c)\+ti,\*(p,c) - tf(p,c)|)j(t) < 0. 
We conclude that 



d 



dt 

Therefore, 



0-£)(*)||li(r") +A*ll(c-c)(t)|| L i (R jv ) ) < C^||(g- ^)(i)|| L i(M^) +^||(c-c)(0IUi(rjv)). 



(e - PXOIU 1 ^) +^ll( c - c)(t)|Ui(]RJV) < e c *(j|£° - q°\\ia(rn) +n\ 



C° - C°||x,i( R JV)j. 



This gives uniqueness, and choosing (f?, c) = (f? m ,c m ), and (g,c) = (0,cb), we find the 
uniform estimates 



(4.6) 



|£m(£)|| L 1 (R iv ); 1 1 Cm 
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As in the purely fluid mechanical model, from the L°° bounds and (4.6), we conclude that, 
for m > 2, 

lbm(t)||/-im^ < Ke ct 



L 1 bounds ON the derivatives OF g m and c m . We differentiate the two equations 
of the system with respect to time, multiply the first one by sign(<9i£> m ) and the second by 
sign((9 i c m ) and use Kato's inequality to obtain 
(4.7) 

d t \d t Q m \ - A(mg™ 1 \d t g m \) < \d t g m \® + g m d p $\d t p m \ + g m d c §d t c m sign(d t g m ), 

d t \d t c m \ - A(\d t c m \) < -d t Q m V sign(<9 t c m ) - g m d p ^d t p m sign(d t c m ) - g m d c ^\d t c m \. 

We integrate in space and add the two equations to obtain, using the monotonicity properties 
of the growth functions, 



d 
dt 



< 



\dtQm(t)\($ + fl^)(0,C B ) + \\Q m \\L«>(Q)\\dc$\\L°°(A)\dtC m (t)\) 



and thus, thanks to assumption (4.1) 



\\9tQ m (t)\\ L i(^if) + \\dtCm(t)\\Li(RV) < Ce (\\(d t g) \\ l i (j& n ) + \(d t c) \\ = Ke 
We can estimate the space derivatives in the same way, and arrive to 



\\dx t gm{t)\\ L l {W N ) + Wd^Cm^WLi^N) < Ke 



Ct 



i = l,...,N. 



Estimates on the derivatives of p m . From the first equation in (4.7), using the strict 



monotonicity of $ with respect to p and the BV estimates that we have already proved, we 
get 



/ \d t g m {T)\ + a If g m \d t p m \ < Ke CT , 
Jr n JJqt 



where, as above, a = min^ \d p $\ > 0. Therefore, 



II 9 t p 



'm\\L 1 {Q T ) 



< 



mg2~ \d t g m \ + 2 



gm \dtPm\ < C(T). 



'Q T n{£i<i/2} J JQ T r\{ Sm >\/2} 

The control on the space derivatives follows by a similar argument, and we obtain 

\\d Xi p m \W{Q T ) < C(T), i = l,...,N. 



Convergence and identification of the limit. The above estimates give strong 
convergence in L\ oc {Q) for g m , cb — c m and p m . Our uniform control of the L 1 norms of 
these quantities in m and t e [0,T] shows that the mass in thin initial strips is uniformly in 
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m small. Therefore, to get global convergence in we just need to control the tails. 

Proceeding as in the case without nutrients, we get 



d 
dt 



Qm(t)<PR < OR' 2 + $(0, C B ) / Qm(t) VR 



from where the tail control for g m , and then for p m , follows easily. 
As for the nutrients, we have 

if (c B - c m (t))(p R < i2 -2 ||A^|| L oo( R N) / (c B — c m (t)) + ^(0,c B )f Qm(t)<PR. 

at Jrn J r n J r n 

But we already have a tail control for g m . Hence, for all R and m large enough, 

/ (c B -c m (t))<pR< f (c B -c°(t))<pR+ f \c°(t)-c m (t)\ip R + CR- 2 + ^(0,c B )e, 
Jm N Jr n Jw n 

which immediately implies a tail control for c B — c m . 
The strong convergence in L 1 (Qt) is enough to pass to the limit and recover the system 



(4.2) in a distributional sense with the limiting pressure graph relation. 



L bounds FOR Vp m , Vc m . The bound for the gradients of the pressures follows from 



equation (4.4) written in the form 

m — 1 



dtPm = "" 2 ^ ^Pm - (TO - 2)|Vp m | 2 + (TO - l)p m $(p m , Cm). 



Indeed, integrating in Q T we obtain 



7TI — J. i 

||Vp m ||i 2(QT) < — — -$(o,c B )||p m || L i ( Q T) + — — -\\p° m \\mR»)- 



This is enough to prove that the limit satisfies the transport equation (4.3). 

The estimate for the gradients of the nutrients, which is needed to prove that the limit falls 
within the uniqueness class, is even easier. We just have to multiply the equation for the 
nutrients by c m and integrate by parts in Qt- Since \1/ > 0, we obtain 



|Vc m || L 2(Q T ) < ||C m || L 2( R JV)/V2. 



5 Uniqueness for the limit system with nutrients 

We meet two difficulties. On the one hand, as in the case without nutrients, p is not a 
Lipschitz, single-valued function of g. On the other hand, the density of nutrients is not 
compactly supported. Hence we can not restrict to the case of a bounded domain. To take 
care of the lack of compact support of c, we will use an idea from [7]. 
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Dual system. Let (f?i,Pi,ci) and (Q2,P2, c 2 ) be two solutions with the same initial data 



and fix a final time T. Since they satisfy (4.2) in the sense of distributions, we have 
{(Qi ~ Q2)d t ip + (pi - p 2 )Atp + (gi$(pi, ci) - Q 2 $(P2, c 2 ))tp} = 0, 
{(ci - c 2 )d t C + (ci - c 2 )AC - (Qi^(pu ci) - £?2^(p 2 , c 2 ))C} = 0, 



(5.1) 




for any pair of test functions 6 T^{Q) such that ip(-,T) = ((-,T) = 0. Denoting 
$i = Ci) and = ^(pi, Ci), we can rewrite the above equations as 

{(gi " 02 + Pi ~ p 2 ) {AdtiJ + BAi; + - CB^) + ( Cl - c^Ity} = 0, 

{(ci - c 2 ) (9 t C + AC - EC) - (Qi -Q2+P1- Pa) (A*! - FS) C} = 0, 

/./Qt 

where 

0<A= <1, 0<5 = <1, 

01 - 02 + Pi - P2 01 - 02 + Pi - P2 

$(pi, Cl ) -$(p 2 , Cl ) $(p 2 , Cl ) -$(p 2 ,c 2 ) 

< C = —£2 < fi, < D = £ 2 < z/ 2 , 

Pi - P2 ci - c 2 

< h = g 2 < ^3, < t — —g 2 < v±. 

Ci - c 2 pi - p 2 

Adding these two equations we get 

(5.2) J J ((Q 1 -Q 2 +p 1 -p 2 )A&,() + (ci-c 2 )B(i>,()) =0, 
where 

.A(V>, C) = + + ( A ®i ~CB)i)- [Mix - FS) C, 
B(i/;,Q = d t ( + A(-E( + D<il>. 

Let G, if be any non-negative functions in T>'(Qt)- If the dual system 

Aty,0=AG, B(ij,Q = H, ^(-,T) = 0, a-,T) = 0, 

admits a smooth solution (ip, () with a suitable decay, we may use ip and ( as test functions 
to obtain 

// ((£?!" £?2)G+( Cl -c 2 )#) =0, 
J JQt 

from where uniqueness would follow. Unfortunately, the coefficients in the equations defining 
the dual system are not smooth. Even worse, A an B may vanish: the system is not uniformly 
parabolic, and a delicate regularization procedure is required to fulfill our plan. 



(5.3) 
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(5.4) 



Regularized dual system. Given C7, H as above, let R > be such that the supports 
of G(-,t), H(-,t), Qi(t), Q2(t), t e (0,T) are contained in Br o (0). Given any R > R + 1, we 
introduce the dual system with regularized coefficients, posed in Q Rj t = B R (0) x (0, T), 

d t ^n,R + ^^n,R + (^l,n ~ ~^^\ ^) nR - (^J 1 n - ^ R = 

dt(n,R + ^Cn,R ~ -^nCn,fi + D n 1p n ^ R = H, 

with final and boundary data given by 

^,k(-, T) = ( n;R (; T) = 0, ^ R = Cn,R = in dB R (0) x (0, T), 

where, thanks to the hypotheses on the solutions, the L 2 (Qt) error in the approximation of 
each coefficient is 0(l/n), and 

l/n< A n ,B n < 1, < C n , D n , E n , F n < M, \®i, n \, \^i,n\ < M, 
\\dtC n \\ L i(Q T ), \\dtF n \\ L i(Q T ), ||V$i ; „|| L 2 ( q t ), || V\l>i in || L 2(Q T) < K. 

Our aim is to use i/j n , R and Cn,_R, suitably extended by zero, as test functions. Since these 
functions have non-zero derivatives at the boundary, we find yet another difficulty: the lack 
of smoothness at the boundary of the ball B R (0), which requires an extra regularization. 

Avoiding the lack of smoothness at the boundary of the ball. For < e < 1/2, 
we consider a family of cut-off functions i] £ e C^°(R ) such that < r\ £ < 1, 

r] £ (x) = 1 if |x| < R — 2e, r\ £ = if \x\ > R — e, 
IIV^IU < C/e, HA^IU <C/e 2 . 



Now put ip = rj £ ip n:R , ( = r] £ ( n ,R in the weak formulation (5.2)-(5.3) to get 

5 

Ve((Ql ~ Q2)G + (Ci - C 2 )H) = -J n>RiE - ¥in,R,s ~ F n>RjE , 



i=l 



where 

Jn,R,e = { (Pi - P2) (2 Vf] £ ■ Vlp n , R + lf} n , R Af] £ ) } , 

" JQr,t 

£n,R,e { (Ci - C 2 ) (2Vt] £ ■ VCn,i? + Cn,R^Ve) } , 

J k,R,e = j j -02+ Pi- P2)VeB n (l ~ {M^n,R ~ Cr$n,R + F n ( n , R ) j 

J n,ij, e = / / {(01- O2+P1- P2)Ve(B - B n ) (A^„ 5J R - C n ^ n . R + F n ( n>R )} , 
J J Qr,t 

ll,R,e = ~ ff {{Ol- Q2+Pl-p2)VeB[^ n , R (C ~ C n ) + ( n , R (F ~ F n )}}, 
J J Qr,t 

In,R,e =jj {(Cl - C 2 ) V£ ( n , R [(E - E n ) + (D - D n )}} , 

J n,R )E = // {(*?! - £2 +P1 - P 2 )%^ (($! - $i, n )Vn,fl ~ (*1 ~ *l,n)Cn,ii)} • 
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The term J n ^ £ vanishes, because the supports of p\(t), P2{t), t e (0, T) are contained in 
Br o (0). Hence we can drop it. 

We will now take the limit e — > 0. The only difficult term at this step is lC n ,R,e- We write 



n,R,e\ <C 

JO J R- 



, VCn,_R \Cn,R\ 
Cl - C 2 | h 



\K 

I R-2e<\x\<R 

Since = on dBn(0), we have 

sup |C„,i?| < 2e|VCn,fl|, N = {R - 2e < \x\ < R, < t < T}, 

N 

and, if v is the unit outwards normal to <9£>r(0), then 



lim SUp |VCn,fl| = SUp |VCn,fl| = SUp \d v Cn,R\ 
£ W \ JV / \x\=R,0<t<T \x\=R,0<t<T 



Therefore, since c\ and c 2 are bounded, 



limsup \K n ,R,s\ < CR sup \d u ( n>R \. 

eiO \x\=R,0<t<T 



Limit n ->■ oo. Uniform bounds. Proving that lim^ooP = 0, % = 1, . . . , 5, is easy, 
thanks to the following estimates. 



Lemma 5.1 There are constants C\, C 2 , depending on T , G and H , but not on n and R, 
such that 

H^rMlll L°°(Qfl,T)> IICrM?IU 00 (Qfl,T) ^ ^1) 



£ 2 (Qb,t) 



<c 2 . 



Proof. L°° bounds. An easy computation shows that the functions 

M(t) = max\ip n>R (x,t)\, N(t) = max\C n>R (x,t)\, 
\x\<R fI<-R 

satisfy the differential inequalities 

-M'(t) + a n (t)M(t) < C(l + M(t) + N(t)) + q n (t)a n (t)N(t), 

-N'(t) <C(1 + M(t)), 

k M(T) = 7V(T) = 0, 

with a n (t) = c ^ Bn (x n , t) > 0, q n (t) = ^r{x n ,t), where x n is the point where the maximum 
defining M(t) is achieved. Thanks to the assumptions on d p § and d p ^, there is a constant 
q > such that < q n (t) < q. 
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Let now 

Q(t) = max (M(t),qN(t)). 
As a combination of the equations on M and N, we find that 

(5.5) -Q'{t) <C{1 + Q{t)), Q(T)=0. 

Indeed, assume that max (M(t ),qN(t )) = M(t ). Then, q n {t )N(t ) < M(t ) and the bad 
a n terms cancel. If, on the contrary, max (M(t Q ),qN(t )) = N(t ), then the result follows 
by the harmless equation on N(t). 

From (5.5) we have Q(t) < 1 — e CT_ *, hence the result. 



Estimates on the Laplacian. We multiply the second equation in (5.4) by d t ( n ,R and 
integrate in space and time to obtain 

II^Cn,fl||l,2(Q JiT ) + ^\\^Cn,R(t)\\h(B R (0)) = [ [ { E n(n,R ~ Dn^n,R + H)d t C,n,R 
z JJQr,t 

< C\\dt(n,R\\L 2 (Q R , T )- 

Hence \\d t Cn,R\\L 2 (Q RT ) * s uniformly bounded. 



We now multiply the first equation in (5.4) by Atp njR — C n ip ntR + F n (n,R and integrate in 
space and time. We get 



^l|V^n,fl(i)||i2 (Bjl( o)) + j t j ^( A VVii - Cnlpn,R + F n (n,R 

= - [ + F^ n , R C n , R ) (t) 

Jb r (o) V 1 J 



, T f ( i>l R dtC n 

+ / / \~ lpn,R(n,Rd t F n + F n 1p n , R d t ( n ,R 



Jt Jb r (o) v j 

Jt JB R (0) V 
+ / / Un,R^G - C n GlP n>R + F n G( n)R ) 

Jt Jb r (o) v j 
<K(l-t + £ 1^^(5)11^(^(0)) dsj . 

Using the Gronwall lemma we get sup t6 ( 0T \ || VVvr(£)IU 2 (.Bk(o)) < C, and then the desired 
estimate □ 

In order to pass to the limit in the /C term we need a good enough estimate, uniform in n, 
for the decay of the normal derivative. 
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Lemma 5.2 As R — > oo, 



(5.6) sup \d„<; n , R \ = oiR 1 -") 

\x\=R,0<t<T 

uniformly in n. 

Proof. As a first step we obtain a more precise estimate for the size of ( n ,R- We first notice 
that, since E n ,D n > 0, D n and H are bounded and compactly supported in B Ro (0), and 
ipn^n is bounded, then 

d t (n,R + AC„,/j = E n C n ,R - D n i) nyR + H> -C1 Brq (0) 

for some C > 0. On the other hand, given a fixed r > 0, if K > is large enough, the 
function 



Z(x,t) = Ke 



-At 



C~ 



|x|V(4(T-t+r)) 



(4tt(T - i + r))^/ 2 
satisfies 



e -\x\y(4(T-t+r)) e -^/(4r) 



Z. + AZ = -Ai^e- (4?r(T _ t + r))JV/2 < - 1 ^o(o)^ (47r(T + r))J v/ 2 
< -Cl Bjlo (o)- 

Moreover, Z(x,T) = 0, and Z(x,t) > for |x| = R and < t < T. Therefore, comparison 
yields 

< ( n ,R < Z in Q RjT . 
If N > 3, we consider a function g = g(x) defined by 

/ \ d 
9( x ) = TZin^ + e > 



\x\ 



where e, d satisfy 



d e -rf»/(4(T+r)) d 

+ e = K— — NAr/a , — — + e = 0. 



(i? - l)^- 2 (4vrr) iV / 2 ' i^" 2 

Note that Ag = on {R - 1 < |x| < i?}. Moreover, 

M =R-i, 0<t <T, 

9{x) = (u,r( x , 0=0, |x| = i2, < t < T, 

g{x) > (u,r( x , T) = 0, R-K\x\<R. 

Therefore, g(x) > ( njR (x,t) for all R - 1 < \x\ < R, < t < T. Since g(x) = = ( n ,R(x,t) 
for \x\ = R, < t < T, we conclude that 

- Cn,*)0M) for \x\ = R, < t < T, 
27 



and hence 

sup \d u (n iR \ < sup \d v g\ 

\x\=R,0<t<T \x\=R,0<t<T 

The result then follows from the following estimate, 
R N -i\d u g\(x) = (N-2)d= [ — 2)KG 



( 47rr )JV/2 \(R- 1)^-2 R N ~ 2 

= o(l) as R —> 00. 

An analogous computation with g(x) = dlog \x\ + e if N = 2, and g{x) = d\x\ + e if N = 1 



lead to the same estimate (5.6). □ 
Limit R -)• 00 and conclusion. Passing to the limit in R, we finally obtain that 

((£?i- £? 2 )G + ( Cl -c 2 )#) =0. 



Hence Q\ = g 2 and c\ — c 2 . Uniqueness for p then follows just taking ip = p 1 — p 2 as test 



function in (5.1 ) 



Appendix: Examples for the purely mechanical model 
A.l Tumor spheroids 

A typical application of the Hele-Shaw equations is to describe tumor spheroids [HI [121 H31 
2U[26l[36]. When nutrients are ignored, the tumor is assumed to fill a ball centered at 0, 

:= { Poo (t) > 0} = { 6oo (t) = 1} = B m (0). 



The radius R(t) of this ball is computed according to the geometric motion rules (2.9) 



and (2.10); that is, we consider the unique (and thus radially symmetric) solution to 
(A.l) - A Poo (t) = $(poo(0) in B R(t)(fy, Poo{R(t),t) = Q, 

and evolve the radius according to 

(A.2) R'(t) = V=\V Poo (R(t),t)\. 

Then, we consider defined as 

( A -3) Qoo(t) = 1b b(1) (o)- 

This is indeed a correct solution to our model. 
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1.0 



0.8- 
0.7- 
0.6- 
0.5- 
0.4- 
0.3- 
0.2- 
0.1- 

o.o- 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

Figure 3: Traveling wave. A traveling wave solution to the mechanical model (1.1), (1.3) 
in one dimension with m = 40. The upper continuous line is g; the bottom dashed line is p. 
Here pu = .85. 



Theorem A.l Let R(0) = R be given. Problem (A.1)-(A.3) defines a unique dynamic 



R(t), g<x>(t), Poo(t), which turns out to be the unique solution to the Hele-Shaw limit prob- 



lem ( |2.3 ) with initial data g^ 



B R o(0)- 



For long times it approaches a 'traveling wave 7 



solution with a limiting speed independent of the dimension, 



(A.4) 



R'(t) — ► y/2Q(p M ) 



t— >oo 



Qip) 



<&{q)dq. 



The limit profile can also be calculated and is one- dimensional. 

For several more elaborate one dimensional models, it is also possible to compute the 
traveling waves which define the asymptotic shape for large times; see |JT] • 

Proof. Existence. A unique solution to (A.l) can be defined for every fixed R(t), be- 



cause the elliptic problem (A.l) comes with a decreasing nonlinearity. Indeed, shooting 



or iterative methods apply and give a solution with smooth dependency on the parameter 



R(t). Therefore the differential equation (A. 2) can be solved thanks to the Cauchy-Lipschitz 
theorem. 

Equivalence with the Hele-Shaw problem. To show that this geometric method gives the 



solution to (2.3), we compute, in the distributional sense, that 

<9t£oo = R'(t)8dQ(t)- 

Also, because p^ satisfies Ap^ = $(poo) in fi(£) and since p^ vanishes out of £l(t), we have, 
still in the distributional sense in R , with v the unit outward normal to O(t), 



Apoo + g^{p 

c 



8dn(t)dvPoc = -$dn(t)\Vpoo(R(t),t)\, 
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and the last inequality follows from the Dirichlet boundary condition. This means that both 



formulations, geometric motion of the domain or equation (2.3), are indeed equivalent 



Asymptotic behavior. We use the notation p R (\x\/R,t) := Poo(x,t). Now p R is solved on 
the unit ball as 

(A.5) - p" R {r) - -^lp' R ( r ) = R 2 $(p R (r)) for < r < 1, p R (l) = 0, p' R (0) = 0. 

By the comparison principle, p R is increasing in R and p R < pm- Passing to the limit in the 
equation we find that p R — > p M as R — > oo. It is easy to see that p' R {l) = O(R), because 

we can build sub- and supersolutions of the form p(r ) = p^[ II — ^-y with A = CR for 

appropriate values of C. By elliptic regularity we conclude that ||p / r ||l°o(r x (o,t)) is of order 
O(R) and, since p' R < 0, it is bounded in L l (M. x (0,T)). As a consequence, ||p'ijL 2 (Rx(o,T)) 
is of order 0(\/R). 



Next we multiply equation (A.5) by p' R and integrate in (0, 1) to obtain 



(P«(l)) a + (N - 1) 



1 (p' R (r)Y 



dr = R 2 Q(p R (0)) w R 2 Q(p 



Mi 



From the above L 2 estimate on p' R , we conclude that p' R (l) ~ R\/2Q(pm)- Going back to 
the dimensionalized variable, (A. 4) follows. 



One- dimensional profile. We consider moving coordinates, s = r — R(t), pt(s) = Poo(s + 
R(t),t). Then pt(s) is positive if and only if s < 0. Hence, 



Vt{s) + 



N-l 

s - R(t) 



p' t (s) + $(p t (s)) = for s < 0, p' t {0) = -R'{t). 



Since R(t) — > oo as t — > oo, we can pass to the limit uniformly on compact sets — K < s < 
to obtain 

$L(s) + $(P°o(s)) for S <0, p / oo (0) = - v / 2gfei 



□ 



Remark. The asymptotic growth rate (A. 4) is generic for solutions to (2.3), since one can 



always sandwich, after a certain delay, any "standard" initial data between two particular 



profiles of the form (A. 3) which asymptotically travel with the same speed (A. 4). 



A. 2 Examples with strong time discontinuities for the pressure 

Here we produce examples of solutions with a jump discontinuity in the pressure as a function 
of time. 

The first and simplest example is constructed in radial symmetry. It consists of an initial 
datum which is constant cq < 1 in a ball centered at and of radius 1 and zero otherwise. 
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Then the solution will have pressure for a certain time, < t < t\, in which the density 



grows exponentially according to equation (2.3), so that 



Qoo(x,t) = c e* (Q) *l Bl (o). 

At t — ti — — Qoo reaches the threshold level g^ — 1 in B x . It cannot happen that 
Poo{t) — for t > t%, since otherwise g^ would violate the bound g^ < 1. Therefore, Poo(t) 
satisfies Ap^ + $(poo) = in some ball B R ^ with R(t) > 0. An easy argument using the 
maximum principle shows that R(t) > 1. It means that Poo{t) is equal or larger than the 
solution poo of the problem Apoo + $(poo) = 0, poo = for |se| = 1, poo > for \x\ < 1. 
Taking the limit t \. ti, we obtain a jump discontinuity in time at t — t\ for every \x\ < 1. 
Note that Poo{t) is a discontinuous function of time with values in any L P (WL N ), p > 1. 

The second example is best presented in one space dimension and consists of the evolution 



of an initial data consisting of two copies of the standard example presented in section A.l 
(tumor spheroid), which have now an interval as support. We locate the supports at a 
distance from each other of say 1. For a time interval, < t < t\ they evolve independently. 
At t — t\, the supports meet and an easy application of the maximum principle, shows that 
the solution is strictly positive in the union of the two intervals, therefore, larger than the 
solution corresponding to the union of the two intervals at t — t\. An easy inspection of the 
solutions at times t± and tf shows that there is a jump discontinuity in the space pressure 
profiles. This example can be adapted to several space dimensions by replacing disjoint 
intervals by disjoint concentric annuli. 



A. 3 The effect of the equation on 

We consider the cell density (defined for short enough times) with two discontinuities 

Qoo (t) = (1 - q(t))t BRi(t) + q(t)t BR2 , R^t) < R 2 , q(t) < 1. 
We claim that the correct dynamics is defined by the speed and values of q(t) given through 

R i® = Y^j\^Poo{Ri(t),t)\, q'(t) = g(t)$(0) while q(t) < 1, 

with poo vanishing outside the ball of radius R±(t) and satisfying — Ap^ = $(poo) in B Rl ^, 
with zero Dirichlet boundary condition. 



These rules are simply derived from equation (2.3). For R\(t) < \x\ < R 2 we have p^ = 
and thus the equation (2.3) is reduced to ^^(t) = ( &(0)^ oo (t) which gives us the evolution 
of q(t). For \x\ < Ri(t), we have d t Qoc = R[(t)(l - q(t))5{\ x \= Rl (t)}, while Ap^ + £oo$(Poo) = 
|Vpoo|&{>]=Ri(t)}- Using the equation (2.3), we get the dynamics for Ri(t). 
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